Quantum Critical Point of Itinerant Antiferromagnet in Heavy Fermion 
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A quantum critical point of the heavy fermion Ce(Rui_ a; Rh a; )2Si2 (x — 0, 0.03) has been studied 
by single-crystalline neutron scattering. By accurately measuring the dynamical susceptibility at 
the antiferromagnetic wave vector = 0.35c*, we have shown that the energy width r(fes), i.e., 
inverse correlation time, depends on temperature as r(fcs) = ci + ciT z ^ ±aA , where ci and C2 are x 
dependent constants, in a low temperature range. This critical exponent 3/2 ± 0.1 proves that the 
quantum critical point is controlled by that of the itinerant antiferromagnet. 

PACS numbers: 71.27.+a, 71.10.Hf, 75.30.Mb, 75.40.Gb 



Quantum critical points (QCP) separating ferromag- 
netic or antiferromagnetic states from paramagnetic 
Fermi liquid states in strongly correlated electron sys- 
tems have been investigated for decades. Successful de- 
scriptions of critical behavior close to QCPs were tradi- 
tionally provided by the self-consistent renormalization 
(SCR) theory of spin fluctuations 0,0] for d-electron sys- 
tems based on the Hubbard model. The mean- field- type 
approximations used in the SCR theory were justified 
by the renormalization group studies [3j based on the 
Hertz effective action above upper critical dimensions. 
For the ferromagnetic QCP, theoretical predictions were 
supported by experimental studies of d-electron systems 
P, . In contrast there is little experimental understand- 
ing of the antiferromagnetic QCP ■ 

A recent intriguing issue of QCP under controversial 
debate is directed toward revealing relevant fixed points 
for antiferromagnetic QCPs in heavy- fermion systems 
In energy scales much lower than Kondo temperature 7k, 
/ and conduction electrons form composite quasiparticles 
with a large mass renormalization in paramagnetic heavy 
fermions. By tuning a certain parameter, e.g., pressure 
or concentration, an antiferromagnetic long range order 
emerges from the Fermi liquid state. In a weak coupling 
picture, it has been hypothesized that the same QCP as 
the d-electron itinerant antiferromagnet, referred to as 
spin density wave (SDW) type QCP, is relevant to the 
heavy fermion QCP 0,0- 

However despite a number of experimental studies of 
heavy-fermion systems showing non-Fermi liquid behav- 
ior, none of them definitely supports the SDW QCP 
0,0 El- This stems partly from experimental difficulty in 
measuring weakly divergent quantities around QCP es- 
pecially for bulk properties, which has been also the case 
for d-electron itinerant antiferromagnets 0] . On the con- 
trary, several recent experiments suggest the possib ility 
of a novel strong coupling picture of the QCP |.4 llfllllj . 
Among these studies direct measurements of the diverg- 



ing spin fluctuation using single-crystalline neutron scat- 
tering for the heavy fermion CeCus.gAuo.i provided in- 
teresting insight [lfjj. On the basis of the observed E/T 
scaling with an anomalous exponent |ld| | and effective 
two space dimensions |l2J. a scenario of a locally criti- 
cal QCP was proposed fl3| . In contrast to the SDW 
QCP, this theory stresses that separation of / spin from 
the quasiparticle state occurs abruptly at the QCP. 

In the present work, we have studied an anti- 
ferromagnetic QCP of another heavy-fermion system 
Ce(Rui_ a; Rh a ;)2Si2 (x = 0, 0.03) using single-crystalline 
neutron scattering. Stoichiometric CeRu2Si2 is an 
archetypal paramagnetic heavy-fermion with enhanced 
C/T ~ 350 mJ/K 2 mol and T K ~ 24 Kill. Extensive 
neutron scattering studies of CeRu2Si2 |l5| have shown 
that spin fluctuations possessing three-dimensional (d = 
3) character are excellently described by the SCR theory 
for heavy fermions 0. A small amount of Rh doping 
x > x c ~ 0.04 induces an antiferromagnetic phase 
[see Fig. ^a)] of the sinusoidally modulated structure 
with the wave vector fc3 = 0.35c* [l7| . 

The Rh doping modifies exchange interactions, while 
keeping Tk constant , and brings about the antiferro- 
magnetic phase in the concentration range x c < x < 0.35. 
The dome structure of this phase and another antiferro- 
magnetic phase with the modulation vector (^0) in the 
higher concentration range 0.6 < x suggests a certain 
frustration effect among exchange interactions. Samples 
nearly tuned to the lowest concentration QCP (x ~ x c ) 
show non-divergent C/T (T — > 0) (lif and Ap cx T 3 / 2 
[see Fig.|Hb)], which are consistent with the SDW QCP 
in d = 3. Thus one can expect that Ce(Rui_ 2; Rh a ;)2Si2 
(x < x c ) is suited to investigate the SDW QCP with- 
out disorder effects. On the other hand, non-Fermi liq- 
uid behavior observed in the higher concentration range 
0.35 < x was reported to be influenced by disorders 
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In order to experimentally show the SDW QCP in 
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d = 3, spin correlation studied by the renormalization 
group theory 0, 0] should be measured directly by neu- 
tron scattering. The theory of the SDW QCP shows that 
the wave-vector dependent susceptibility for the tuned 
sample (x = x c ) diverges as x(fc 3 ) oc T -3 / 2 0, Q, or 
the characteristic energy of the spin fluctuation, i.e., 
the inverse correlation time, depends on temperature as 
r(fc 3 ) oc x(fc 3 ) _1 oc T 3 / 2 . By taking the detuning effect 
(x < x c ) into account, leading two terms of Tjk com- 
puted by the renormalization group theory 0, l2lj are 
given by 



r(fe 3 ) = ci + c 2 T 



3/2 



(1) 



where c\ and C2 are x dependent constants. This equation 
is an approximation in the temperature range Tfl <C 
T <C 7k, where Tfl is a crossover temperature below 
which system shows the Fermi liquid behavior 0, l2l| . 
In the present work, we have accurately measured r(fe 3 ) 
and have shown that it agrees well with Eq. for both 
the nearly tuned sample x = 0.03 and the stoichiometric 
sample x — 0, which indicates that disorder does not 
influence the critical behavior. 

Neutron-scattering measurements were performed on 
the triple-axis spectrometer HER at the Japan Atomic 
Energy Research Institute. It was operated using final 
energies of Ef = 3.1 and 2.4 meV providing energy res- 
olutions of 0.1 and 0.05 meV (full width at half maxi- 
mum), respectively, at elastic positions. Single crystals 
with a total weight of 19 g (x = 0) and 17 g (x = 0.03) 
were grown by the Czochralski method. Two sets of mul- 
ticrystal samples aligned together were mounted in a He 
flow cryostat so as to measure a (hOl) scattering plane. 
All the data shown are converted to the dynamical sus- 
ceptibility. It is scaled to absolute units by comparison 
with the intensity of a standard vanadium sample. We 
note that a new point of the present work is unprece- 
dented experimental accuracy in determining the critical 
exponent [3/2 in Eq. J[J] using large samples and long 
counting time. This has enabled us to determine the sin- 
gularity of QCP and to make qualitative conclusion of 
the universality class. In the pioneerin g w ork using the 
related compound Cei_ 2; La 2 ;Ru2Si2 |8|,|22j, Eq. (JIJ has 
been discussed assuming the fixed value of the exponent, 
which could be determined only with an order larger ex- 
perimental error. 

The imaginary part of the dynamical susceptibility at 
Q = fc 3 + q with small \q\ and \E\ is predicted to be 
approximated by 



Imx(fc 3 + Q, E) 
T(k 3 + q) 



(2a) 



X (fc 3 )r(fc 3 )£ 
E 2 + T(k 3 + qf ' 
D c [ K 2 c + q 2 c +F{q 2 a +q 2 )] , (2b) 



where D c and F are T independent parameters, and 
k c is inverse correlation length along the c axis |23| . 
This expansion form with the T independent product 
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FIG. 1: (color online) (a) The phase diagram and Tk of 
Ce(Rui_ I Rh a; )2Si2 are reproduced from Refs. 0, 0|. (b) 
Resistivity of the sample with x = 0.03 is plotted as a function 
of T 3/2 . Constant- £ scans taken with E = 0.2, 0.4, 0.6 and 
1.2 meV along (c) Q = (H, 0, 0.65) and (d) (1,0, L) lines for 
the sample with x — 0.03. Data of T = 1.5 and 4 K are 
shifted by 0.25 and 0.1 emu/mol Ce, respectively for clarity. 
Curves in (c) and (d) are calculations using Eqs. 12aH and 12 bL 
corrected for resolution functions with the same fit parameter 
as those shown in Fig. |5] 



x(fe 3 )T(fc 3 ) has been used in the SCR 0, El and renor- 
malization group |2*H theories. The two parameters D c 
and F were determined by usin g co nstant-Q and -E scans 
for both samples with x = [15j and 0.03 at T = 1.5 



K. These data were fitted to Eqs. (|2a|l and l|2b|l con- 
volved with the resolution functions. In Figs. He) and 
^d) we show constant-i? scans through the antiferro- 
magnetic wave vector Q — (101) — fc 3 and the fit curves 
for the sample with x = 0.03. The good quality of fitting 
indicates that Eqs. lll'all and l|2b(l well describe the exper- 
imental data at T = 1.5 K. We obtained concentration 
independent values of the parameters D c = 98 ± 4 (meV 
A 2 ) and F = 0.12 ±0.01. 

Temperature dependence of T(fc 3 ) = D c n 2 c has been 
determined by performing constant-Q scans taken at 
Q = (101) — fc 3 with much higher statistical accuracy 
than previous measurements [l5l |22| . These scan data 
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were fitted to Eq. i12a[| convolved with the resolution func- 
tions, where there are two adjustable parameters r(fe 3 ) 
and x(^s)- Several fit results of the constant-Q scans for 
the samples with x = and 0.03 are shown in Fig. [21 
From these figures one can see that the quality of fit- 
ting is excellent. We also checked the T independence of 
the parameters D c and F by comparing the constant-^ 
scans in Fig. at T = 4 and 8 K with those calculated 
using the T dependent r(fe 3 ) and x(fe 3 ) determined by 
the constant-Q scans. The calculated curves in Fig. ^ 
which have no adjustable parameters for T = 4 and 8 
K, agree reasonably well with the observations. Thus we 
conclude that the theoretical approximation of Eqs. iffll 
has been experimentally confirmed, and that the fit pa- 
rameter r(fc 3 ) has been determined very precisely. 

Temperature dependence of r(fc 3 ) is shown in Fig. |3| 
by plotting data as a function of T 3 / 2 . At low tempera- 
tures observed data clearly agree with the linear behavior 
of Eq. (JJJ. In fact, by least squares fitting we obtained 
r(fe 3 ) = (0.67 ± 0.01) + (0.0095 ± 0.0021)T 1 - 53±0 08 (in 
units of meV) inl.5<T<16K for the sample with x = 
0andr(fc 3 ) = (0.129 ±0.007) + (0.020 ±0.003)T 149±007 
in 1.5 < T < 8 K for x = 0.03. Therefore we conclude 
that the observed critical exponent 3/2 ± 0.1 is in agree- 
ment with the theoretical value 3/2, and consequently 
that the temperature dependence of the spin fluctuation 
is controlled by the SDW QCP in d = 3. The same ex- 
ponent for both x = and 0.03 samples ensures that the 
randomness due to Rh doping does not affect the criti- 
cality. The temperature independence of x(fc 3 )r(fc 3 ) was 
also confirmed as shown in the inset of Fig. |3J 

The constant c\ in Eq. 0} is proportional to the the- 
oretical tuning parameter, a coefficient of the quadratic 
terms of the Hertz effective action, and c\ oc x c — x is 
normally assumed |3J • This assumption is consistent with 
the observed values of c\ and the critical concentration 
x c = 0.04 ±0.005. In contrast to this agreement, the con- 
centration dependence of C2, c-2.{x = 0.03)/c2(a; = 0) ~ 2, 
may suggest a certain difficulty in the theoretical in- 
terpretation. The constant C2, being proportional to 
the coefficient of the quartic term of the Hertz effective 
action, is usually supposed to vary slowly in the con- 
centration range of interest 0|. The appreciable varia- 
tion of C2 seems to suggest unknown perturbations for 
Ce(Rui_ x Rh x ) 2 Si2. In terms of the SCR theories, the 
variation of ci may be accounted for by the adjustable 
mode-mode coupling constant (oc oi) that is used in the 
phenomenological SCR theory developed for d-electron 
systems 0, yj- Despite this problem, we think that the 
critical exponent 3/2 determined by basic characteristics 
of the system, the space dimension d — 3 and the dy- 
namical exponent z — 2, is more important and decisive 
to conclude the nature of the QCP. 

An advantage of the present neutron scattering study is 
that Eq. Q holds in a wider temperature range compared 
to those of indirect measurements using bulk properties, 
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FIG. 2: (color online) Constant-Q scans measured at the 
antiferromagnetic wave vector Q — (101) — &3 for the samples 
with x = (a) and x — 0.03 (b). Curves are fit results using 
Eq. I2al with two adjustable parameters E(fca) and x(^s)- 

e.g., C/T = 70 - oT 1 ' 2 or Ap oc T 3 / 2 H [3. Theo- 
retically Eq. Q is an approximation in the temperature 
range Tfl <§C T -C 7k, where Tfl is the crossover tem- 
perature to the Fermi liquid state 0, . The tempera- 
ture range in which Eq. is observable can be discussed 
quantitatively using the SCR theory [|| . In Fig. the 
dashed line reproduces the SCR computation of r(fe 3 ) 
for CeRu2Si2 based on the previous neutron scattering 
study [15j. Apart from discrepancy of the coefficient C2, 
one can see that the T 3 / 2 dependence of Eq. (JTJ is a 
good approximation for the SCR curve in the T range 
2.5 < T < 13 K (4 < T 3 / 2 < 47), which agrees with that 
of the observed data for CeRu2Si2- Since the lower bound 
temperature is shown to be proportional to the tu ning 
parameter by the renormalization group theory 0, l2l| . 
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FIG. 3: (color online) Energy width r(fcs) of the Lorentzian 
form Eq. 12aH is plotted as a function of T 3/2 . Full lines arc 
fit to r(fc;i) = ci + ciT v with adjustable parameters ci, cz, 
and v in low temperature ranges, where data are displayed by 
open symbols. The dashed line is the calculation [15| using 
the SCR theory [g. The inset shows temperature dependence 
of the product x(k3)T(ks). 



the T 3 / 2 dependence can be expected in a T range of 0.5 
K < T (0.35 < T 3 / 2 ) for the sample with x = 0.03. The 
smaller T range of the T 3 / 2 dependence for x = 0.03 is 
probably related to the larger constant c^ix — 0.03). We 
note that below 2.5 K the SCR computation of T(ks) for 
x = 15j is approximated by the Fermi liquid behavior 
of r(Jfc 3 , T) - r(fe 3 , 0) oc T 2 [13, which is not clearly seen 
within the present experimental error. 

In connection with neutron scattering experiments of 
CeCus.gAuo.i 0,0, it was theoretically predicted 
that the locally critical QCP is relevant for the two- 
dimensional spin fluctuation, in agreement with the ex- 
periments of CeCu5.9Auo.1- This theory also predicted 
that the SDW QCP is relevant for the three-dimensional 
spin fluctuation, which is in accord with the present re- 
sults. Finally we note that the present work is first clear 
experimental verification of the SDW QCP to our knowl- 
edge among single-crystalline neutron scattering studies 
on QCP or non-Fermi liquid behavior of heavy fermions, 
e.g., Refs. EElEililEiEil and d-electron systems, 
e.g., Refs. |23,|22j. Assuming that criticalities of QCPs 
are classified into a limited number of universality classes, 
we expect that the SDW QCP remains to be observed in 
other systems. 

In summary, we have demonstrated that the quantum 
critical behavior of the heavy fermion Ce(Rui_ 2; Rh a ;)2Si2 
is controlled by the SDW type QCP in three space dimen- 



sions. The inverse correlation time, i.e., energy width 
r(fc 3 ) of the dynamical susceptibility, shows the T 3 / 2 
dependence predicted by the renormalization group and 
SCR theories. 

We wish to acknowledge B. Fak, J. Flouquet, T. 
Taniguchi, and T. Moriya for valuable discussions. 
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